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Abstract 

Recently the charm meson was observed to have a clear branching ratio into 
the low energy vr — vr sigma resonance, while this channel was not detected in the 
decay. It is shown that this is consistent with the standard treatment of exclusive 
charm meson decays and a proposed glueball/sigma picture. 
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The Fermilab E791 Collaboration has recently reported seeing a clear signal of the decay 
of the into a vr"^ plus a low-energy scalar-isoscalar vr — vr resonance, with a (mass,width) 
of about (480,325) MeV[0. This is presumably the elastic resonance found]^ in the analysis 
of TT — TT scattering with (mass, width) of about (400,400) MeV. We refer to this as the sigma, 
cr. In the D+ decays, however, the sigma was not found0. It is the purpose of the present 
work to show that these experimental results are consistent with the glueball/sigma model 
which we have recently proposed. With this model there is a unique skeletal diagram which 
allows us to calculate the ratio of the (T7r+ to the 07r+ decay widths of the D"*" with no new 
parameters. The QCD corrections to the skeletal diagrams for these decay channels are then 
discussed. For the decay into the (J7r+ channel it is shown that there are two skeletal 
diagrams then tend to cancel, which can explain why this channel has not been observed in 
experiments. 

The sigma glueball idea is based in part on the QCD sum rule analysis of the scalar 
gluon correlator and in part from the large a decay rates of scalar glueball candidates. 
Although the QCD sum rule method is not used in the present work, we briefly review 
the method for finding glueball masses. Using the composite field operator for the scalar 



glueball, J'^(x) = a^G^, where is the gluon field tensor and J'^{x) is proportional to the 
pure-glue term in the QCD Lagrangian, the scalar gluon correlator is defined as 

n(p) = i j d^x e'^" < I T[J{x)J{Q)] I > , (1) 

where the link operator has been omitted. The correlator is evaluated numerically in lattice 
gauge calculations. In the QCD sum rule method the dispersion relation for the correlator 
is equated to a QCD operator product expansion, with the mass of the Glueball determined 
by estimating the value of the pole in the dispersion relation. Using a subtracted dispersion 
relation and carrying out the standard QCD sum rule analysis[0] a number of theorists 
find a glueball solution in the region of 300-600 MeV|^, depending on the values used 
for the nonperturbative condensate parameters. In a recent study of the coupled scalar 
mesons and glueballs in which instanton as well as other gluonic effects have been included P| 
we find a mainly meson solution corresponding to the fo(1370), a mainly glueball solution 
corresponding to the fo(1500), and a light glueball in the region of 400-500 MeV^], consistent 
with our earlier studies. The main two sources of error in the method are the determination 
of the values of the condensates and the treatment of the continuum part of the dispersion 
relation. A light glueball has not been found in quenched lattice gauge calculations. 

The other observations which lead to the glueball/sigma picture are first that the fo(1500), 
which has characteristics of a glueball with a scalar meson admixture, has the four-pion 
channel as its largest decay branching fraction|]^; and second, that a BES analysisj^ has 
shown that the four-vr channel is dominated by two-sigmas. This is also true for the fo(1710) 
and fo(2100) which have recently been shown to have glueball characteristics. This has led 
to our conjecture that the light glueball and sigma form a coupled-channel system, with 
the sigma resonance driven by the glueball pole. From the Breit-Wigner resonance fit to 
the resonance 0] we obtain the matrix element < a \ T^*"* | GB > for the coupling of the 
scalar glueball to the scalar tt — tt resonance. Note that the coupling interaction is not 
needed in this work, only the matrix element taken from the experimental analysis. We 
have used this for the study of hybrid baryon decay [^], and the pomeron-nucleon coupling 



and the production of sigmas in high energy p-p collisions via the pomeron||10||. This is the 
glueball/sigma model. Since it is possible that the treatment of the gluonic continuum could 
give a false solution for the light glueball, we consider this picture to be a conjecture, however, 
the model would also follow from a very broad gluonic structure in the region of the a (i.e., 
the glueball pole is far from the real axis) which couples strongly to the sigma resonance. If 
the picture is valid there are many important consequences that can be observed in a variety 
of experiments. We believe that the charm meson decays are an excellent example, which is 
the subject of the present work. 



It is convenient use the quark-diagram classification schemelTTI], in which there are six 



skeltal diagrams, by which we mean processes without explicit gluonic effects. The need to 
consider processes in addition to what was once considered to be the dominant spectator 
decay process was made evident by the large difference between D"*" and D° lifetimes 



See Ref [|T^ for a review of exclusive D"*", D° and D+ charm decays expressed in terms 



of this classification. Theoretical treatments have been mainly based on an effective weak 
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Fig. 1 decay to tc"*" and (|) or o via internal W"*" emis son 



Hamiltonian based on the standard model (see Refs.|T^, |T5l for reviews of the method). For 
example, for the decays of the D"*" to a + vr"*" or a + 7r+, which are of central interest for 
the present work, the effective Hamiltonian is 



(2) 



where {qe(lf)L = qtJ^lq'f = Qel^i^ ~ 1^)1% ^ij the CKM matrix elements and the 
color is summed. The constants Ci,C2 have been estimated by renormalization group 
calculations ||16[. Calculations with this framework proved to be quite successful for the 



study of most exclusive D decays. During the period when the method was being devel- 
oped, and many exclusive charm decays were being measured, however, it was argued that 
additional flavor singlet gluonic processes, called hairpins, might be important ill], O. It 



was also shown that final state interactions must be considered on the same level as the 
hairpins|T8|, and that strong interaction effects at least for certain exclusive decays make it 
difficult if not impossible to identify the contributions of the varous skeletal quark diagrams. 
For inclusive processes one can carry out operator product expansions and treat the non- 



perturbative QCD effects using known condensates (see, e.g., Ref [0] for a recent review), 
however, for exclusive processes accurate inclusion of QCD is difficult. 

For the decays D+ ^ 07r+ and D+ ^ o"7r+ there is a unique skeletal diagram, the 
internal mechanism illustrated in Fig. 1. Using Eq.(^ with a Fierz transformation of 
the second term one has 



< X7r+ I H^^^ I > 



—j=cos6cSin6cC-i < X | ss)l | >< tt" 



{uc)l \D+> . (3) 



The constant C3 will not be needed here. The matrix elements of the quarks currents for 
X = (h^a are 



< I sJ^s 



< I sJ^s I a > 



(4) 
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with /o- the couphngs given by the short-distance values of the wave functions of the 0, cr, 
respectively. For the decay — > crvr^ there is additional contribution from the process 
with the s and s from the two W+ vertices being replaced by d and d. The CKM matrix 
elements are about the same and < | dJ^d \ a >~ 1.2 < | sJ^^s | a >, since the s 
quark condensate is about 80% of the d quark condensate. Thus the effective value of fa is 
increased by about a factor of two. 

Here we use for the value of the cj) coupling[|l^ = 228 Mev. The glueball amplitude 
constant can be estimated from the low-energy theorem^ 

n(0) = - < I I > , (5) 

with < I I > the gluon condensate. The cr coupling can be obtained from another low 
energy theorem 0] 

i j d^x<Q \ T[ss(x)a,G2(0)] | > ~ 247r < | | > /6o , (6) 

with ho = 9.667 for three colors and two flavors, and Eq.(^, with a similar form for the 
d-quark contribution to give 

< I sJ^s I (J > + < I dJf,d \ a> = O.USf,oGeV^. (7) 

Taking the sigma mass to be 500 Mev, the phase space ratio without spin a/(j) = 1.23. 
Therefore from Eq.(|^) we obtain for the ratio of the exclusive decay widths 



r(D+ ^ an 
T{D+ (f)n 



- - 0.15. (8) 



The experimental value of this ratio from Ref.[jl| for the a and the PDG0| for the 
is 0.22. In the light of both the uncertainties in the calculation and the possibilitiy of 
other contributions, discussed next, there is quite reasonable agreement between theory and 
experiment. 

Although the diagram shown in Fig.l is the only skeletal diagram that contributes to 
either the or a decays of the D"*", there are two other diagram with explicit QCD processes 
that also contribute: the annihilation-hairpin processes (Fig. 2a) and the spectator-final state 
interaction process (Fig. 2b). The later is analogous to the process treated in Ref. [0 in 
which it was shown that such final state rearrangement processes might be comparable in 
magnitude to the nonspectator decays. We do not attempt to estimate these prosesses which 
are very difficult to calculate accurately. 

For the decay into a 07r+ there is one skeletal diagram, the spectator shown in Fig. 3a; 
while for the decay into a crvr"^ there is a contribution from the spectator graph and also two 
annihilation graph processes, shown in Fig.3b. For the — > crvr"'' the matrix element of 
the weak effective Hamiltonian, after a Fierz transformation, is 

< an+ I H^ff I > = ^cosOliCi + 3C2) < (Jn+ \ {ud)L | >< | {sc)l \ > (9) 

v2 
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Fig. 2 (a) Annihilaton-hairpin, (b) Final state interacton 
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Fig. 3 (a) Spectator, (b) Annihilation 
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Fig. 4 Diagram for T — > 71 a (or scalar meson) 



Using the results given in Ref. one sees that the ratio of the a to (p rates for the 
decays, which contains a factor of (Ci + 3C2)^/C^, would be strongly reduced compared to 
the D+ decay. This gives a qualitative explanation for the fact that the o"7r+ channel was 
not found 1^ in the E791 experiment. 

It is interesting to note that the only skeletal diagram for the + 0(cr) decays is 

the annihilation process; and that there are also two D° decays to the a or 0, D" (f){a)K° 
and D° — (f){a)K°, that have the W-exchange as the only skeletal diagram. We expect that 
the ratio of the a to the fraction should be similar to the corresponding D"*" decays described 



above. Another interesting related experimental observation by the CLEO Collaboration! 21 
is that in the r~ — >■ z/^7r~7r°'7r° decay the data cannot be fit without a sigma channel. The 
diagram with the W~ coupling to the an~ through a du, illustrated in Fig. 4, leads to this 
channel in a natural way in the glueball/sigma model. In this case a more detailed theoretical 
study of the {qaQb) l matrix elements is needed to obtain the branching ratios, which we do 
not attempt here. 

We conclude that the D"'' and decays are consistent with our glueball/sigma model 
and suggest that experimental searches for the sigma channels in other charm decays, as well 
as other heavy quark and lepton hadronic decays, would be rewarding. 

^ On leave from Carnegie Mellon University 
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